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We describe a search for a third generation vector leptoquark (VLQ3) that decays to a & quark 
and tau lepton using the CDF II detector and 322 pb -1 of integrated luminosity from the Fermilab 
Tevatron. Vector leptoquarks have been proposed in many extensions of the standard model (SM). 
Observing a number of events in agreement with SM expectations, assuming Yang-Mills (minimal) 
couplings, we obtain the most stringent upper limit on the VLQ3 pair production cross section of 
344 fb (493 fb) and lower limit on the VLQ3 mass of 317GeV/c 2 (251GeV/c 2 ) at 95% C.L. 

PACS numbers: 14.80.-j, 13.85.Rm 
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Despite its extraordinary success, the standard model 
(SM) of elementary particles has structural deficiencies. 
The parallels between the families of quarks and lep- 
tons suggest a possible link between these two sectors 
at higher mass scales. Leptoquarks, therefore, have been 
proposed as fractionally-charged color-triplet bosons car- 
rying both lepton and baryon quantum numbers. Lep- 
toquarks appear in a wide range of theories, including 
SU(5) grand unification [H, superstrings 0, SU(4) Pati- 
Salam [3j, and compositeness models |4j; direct searches 
for leptoquarks at hadron colliders probe the parameter 
space of some of these models. 

The various leptoquark states are classified according 
to the quantum numbers of SM gauge group interac- 
tions At the Tevatron collider, these states would be 
predominately pair produced through quark anti-quark 
annihilation. In general, larger cross sections, and thus 
better search sensitivities, are predicted for vector (spin 
1) than for scalar leptoquarks |6j. 

We search for third generation vector leptoquark 
(VLQ3) pair production, and assume each VLQ3 decays 
promptly to a b quark and a tau lepton. As the tri- 
linear and quartic couplings between vector leptoquarks 
and gluons can have model-dependent "anomalous" con- 
tributions, we examine two scenarios: one with Yang- 
Mills couplings, where vector leptoquarks are fundamen- 
tal gauge bosons of an extended gauge group, and the 
other with minimal anomalous couplings p, qJ . Previous 
VLQ3 searches have been carried out in experiments at 
pp, e + e~, and ep colliders [1, 0, [1(1 11 1. Our new re- 
sults substantially extend the reach beyond the previous 
limits. 

The results reported in this Letter are obtained from 
data corresponding to 322 pb -1 of integrated luminos- 
ity collected between March 2002 and August 2004 by 
the CDF II detector [l2j], operating at the Tevatron pp 
collider. Detector components relevant to this search 
are described briefly here. The charged particle track- 
ing volume, which is inside a uniform 1.4 T magnetic 
field, includes a multi-layer silicon microstrip detector 
system and an open-cell drift chamber (COT) that pro- 
vide position, momentum, and charge information in the 
range \r)\ < 1.0 [l3j]. Outside the solenoid, electromag- 
netic (EM) and hadronic calorimeters are arranged in a 
projective tower geometry to measure electron, photon, 
and jet energies. A set of strip and wire chambers (CES) 
is embedded in the EM calorimeter, at the depth where 
the longitudinal development of EM showers is expected 
to be maximal, and provides information used for the 
identification of electron candidates and reconstruction 
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of 7T° candidates that decay through ir° — > 77. Outside 
the calorimeters, drift chambers and scintillators provide 
muon candidate identification in the range \rj\ < 1.0. 

This search assumes a branching ratio B(VLQ3 — > 
br) = 1 , and considers a signature where the decay prod- 
ucts of the VLQ3 pair, T + T~bb, yield two jets from the b 
quarks, an electron or muon from a lcptonically decaying 
tau, and a hadronically decaying tau (r^). We do not 
attempt to identify the jets as originating from b quarks, 
as we find this would degrade the search sensitivity. A 
three-level trigger system selects events with lepton can- 
didates and charged tracks [l4(. These events are then 
classified as er^ or fiTh based on the flavor of the leptonic 
tau decay. 

Selected events are required to contain at least one 
well-identified electron (muon) candidate that passes 
fully through the fiducial volume of the COT, with 
transverse energy (momentum) Et > lOGeV (pr > 
lOGeV/c) [l3j| . To reduce the background due to multi- 
jet quantum chromodynamics (QCD) events, an isolation 
requirement is imposed upon the electron or muon can- 
didate. Specifically, the sum of the pr of all additional 
tracks within a cone in -q — <p space of AR < 0.4 around 
the track direction of the candidate is required to be less 
than 2 GeV/c, where AR = ^ ' Aff + A0 2 . 

Hadronic tau candidates are formed by matching nar- 
row clusters of calorimeter towers with tracks. The proce- 
dure is described briefly here and in detail elsewhere [l5[ . 
A cluster is seeded by a tower with Et > 6GeV. The 
highest px track with px > 6 GeV/ c that points to the 
cluster is defined as the seed track. Other tracks are 
added if they are within an angle of 0.05 to 0.17 radi- 
ans from the seed track, where the angle depends on the 
total cluster energy. Tau candidates with one or three 
tracks are considered. We reconstruct 7r° candidates as 
single narrow strip and wire clusters in the CES, and in- 
clude them as part of the candidate if they are within 
0.17 radians of the seed track and have Et > 1 GeV. 
The Th candidates are required to have \r)\ < 1.0 and 
Et > 15 GeV. To reduce instances where a jet is misiden- 
tified as a 77, , we place requirements on the mass formed 
by constituent tracks, m tr k < 1.8GeV/c 2 , and the mass 
formed from tracks and ir° mesons within the tau can- 
didate, wtrk+TT < 2.5 GeV/ c 2 . To better discriminate 
hadronic taus from electrons, we require the ratio of Et 
deposited in the hadronic calorimeter to the sum of the 
magnitudes of the transverse momenta of the tau tracks 
to be greater than 0.1. To obtain isolated Th candidates, 
we require that the region outside the track selection 
cone, but inside a cone of 0.52 radians from the seed 
track, contains no tracks with pt > 1 GeV/c. In addi- 
tion, the region outside the ir° inclusion cone, but inside 
a cone of 0.52 radians from the seed track, is required 
to contain less than 0.6 GeV of summed Et due to ir° 
candidates. 

To ensure efficient event reconstruction, the electron 
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or muon candidate (I) direction must be separated from 
the tau candidate direction by AR(rh,i) > 0.7. Jet can- 
didates, with Et > 15 GeV, are identified in the region 
\r]\ < 2.4 and are required to be separated from the lepton 
candidates by AR(£/T h ,jet) > 0.8. 

There are a number of SM processes which can mimic 
the VLQ3 signal. The first category consists of back- 
ground processes which contain a real er/, or jJTh plus 
two jets. The primary processes are Z°/j* — > tt plus 
two jets, and tt — > WbWb, where one W yields a hadronic 
tau via W — > r^ T , the other 14 7 similarly yields an elec- 
tron (muon) or leptonically decaying tau, and the two b 
quarks give jets. 

The second category of backgrounds consists of those 
that include misidentified final state particles. These 
include tt — > WbWb, where a jet from a hadronic W 
decay can be misidentified as an electron (muon). The 
processes it — > WbWb, Z /^* — > e + e~ plus jets, and 
Z°/-f* — > plus jets also contribute as backgrounds 

when an electron (muon) or jet is misidentified as the t^. 
Events with VP plus jets can pass the selection if one of 
at least three jets is misidentified as the 7V Contribu- 
tions from diboson channels (WW, WZ, and ZZ) plus 
jets are negligible. The above contributions and their un- 
certainties are estimated using PYTHIA [l6[ Monte Carlo 
simulation and GEANT 17;] CDF II detector simulation. 
Background from multi-jet QCD can contribute when 
jets from quarks are misidentified as an electron (muon) 
or Th- Photon plus jets background enters when high- 
er photons convert within the detector and at least one 
of the resulting electrons appears as a primary electron 
candidate, while a jet is misidentified as the Th- Contri- 
butions from both of these sources are estimated directly 
from the data, using methods described elsewhere 15 1. 

Further event selection reduces the backgrounds. 
Backgrounds associated with misidentification are re- 
duced by approximately one-half through a requirement 
of oppositely charged electron (muon) and Th candidates, 
where the charge of the Th candidate is defined as the sum 
of the charges of all constituent tracks. Events consistent 
with photon conversions and cosmic rays are removed. To 
reduce contributions from Z boson production, events are 
rejected if 76 < m(£,X) < 106GeV/c 2 , where X can be 
the tau candidate or a second electron candidate in the 
erj, channel, or a second muon candidate in the /iTh chan- 
nel. We require fr > 10 GeV, where fr is the magnitude 
of the missing transverse energy fl3} | . This requirement 
reduces backgrounds from Z° /-f* — > r + r _ and multi-jet 
QCD processes and is nearly 100% efficient for the sig- 
nal process. We define Ht as the scalar sum of electron 
(muon) candidate Et, tau candidate Et, event pr, and 
the transverse energies of the two highest Et jet can- 
didates. The requirements on Ht are given below. The 
final selection requirement is that the event must contain 
two or more jets. 



To simulate VLQ3 pair production and decay 18j, we 



have added the production and decay processes to the 
GRACE [l9| matrix element event generator, which cal- 
culates amplitudes, and to the GR@PPA [20j interface, 
which speeds up computations of the interactions of the 
primary hadrons. In addition to providing the theoret- 
ical cross section, these programs yield events that are 
processed through tauola 2l| to simulate tau decays, 
PYTHIA [l(| to simulate parton showering, fra gme ntation, 
and additional particle decays, and GEANT [17J] for the 
full CDF II detector simulation. For the first time, this 
framework includes helicity amplitudes for the full matrix 
element at tree level and propagation of helicity informa- 
tion from the leptoquarks to the tau decay products. We 
use the parton distribution functions (PDFs) CTEQ5L [22| 
and renormalization energy scale Q 2 = ttivlqs 2 - 

We determine total selection efficiencies, including fac- 
tors for triggering, geometrical and kinematic acceptance, 
particle candidate identification and isolation, and back- 
ground suppression criteria. The total efficiencies, av- 
eraged between the er/, and /j,Th channels, range from 
about 2.2% (1.4%) to about 6.1% (5.9%) for Yang-Mills 
(minimal) couplings over the mass range 160GeV/c 2 to 
400 GeV/c 2 . For the example ofm VLQ3 = 320 GeV/c 2 and 
Yang-Mills couplings, the efficiencies for the eTh and /iTh 
channels are (6.0 ± 0.1)% and (6.1 ± 0.1)% respectively. 

We define two signal regions. In addition, three side- 
band control regions in the plane of the number of jets 
(TVjets) versus Ht are used to verify the expected compo- 
sition of the backgrounds and the distributions for kine- 
matic quantities. The primary signal region (SRa) has 
^Yjcts > 2 and Ht > 400 GeV, and is sensitive to the 
highest mass leptoquarks. The secondary signal region 
(SR B ) has N icts > 2 and 250 < H T < 400 GeV, and 
adds sensitivity to lower VLQ3 masses (down to the pre- 
vious mass limits). The three control regions are called 
CR0J, CR1J, and CR2J, where 0J, 1J, or 2J specifies 
the number of jets (0, 1, or > 2). Regions CR0J and 
CR1J include the Ht range Ht > 80 GeV, while region 
CR2J is restricted to 80 < H T < 250 GeV. Table U shows 
the expected background contributions in the control and 
signal regions, as well as the number of events observed 
in all regions. The signal regions are examined only after 
an a priori optimization of the Ht ranges that maximizes 
signal sensitivity. Control region CR1J in the [iTh chan- 
nel contains the largest difference between the expected 
and observed number of events, with a difference of 1.9 
sigma. Figure Q] shows the Ht distributions for the er/j 
and pLTh channels, and includes the control region CR2J 
and the two signal regions. 

The dominant sources of systematic uncertainties on 
the signal efficiencies are the amount of initial state radi- 
ation (ISR) and final state radiation (FSR), the tau iden- 
tification, and the isolation requirements. The ISR and 
FSR uncertainties, as evaluated by varying the amount 
of ISR and FSR in simulation, are each approximately 
3.7% of the selection efficiency. The tau identification 
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TABLE I: Numbers of events observed in data and estimates 
for the total background, for the er^ and fiTh channels, in 
the control regions (CROJ, CR1J, CR2J) and signal regions 
(SRb, SRa). For the backgrounds, the statistical uncertainty 
is given first, followed by the systematic uncertainty. 







er h 




fJ-Th 




Data 


Background 


Data 


Background 


CROJ 


129 


122.1 ±2.1 ± 11.3 


129 


147.1 ±2.6 ±12.3 


CR1J 


110 


109.2 ± 2.3 ± 9.3 


79 


100.5 ±2.5 ±6.7 


CR2J 


36 


33.4 ± 1.4 ±4.8 


26 


30.6 ± 1.6 ±3.8 


SR B 


5 


3.3 ±0.3 ±0.5 


3 


2.2 ±0.3 ±0.3 


SRa 





0.3 ±0.1 ±0.1 





0.2 ±0.1 ±0.1 



f Data (\Ldt=322pb'') 

VLQ3VLQ3 (m=320 GeV/c 2 j 

QCD 




100 200 300 400 500 600 700 800 900 1000 

^(GeV) 



FIG. 1: Distributions of Ht for the (a) eru channel and (b) 
fiTh channel, including one of the control regions used for 
validation (80 < Ht < 250 GeV), and the regions used for 
the search (Ht > 250 GeV). The individual background con- 
tributions are shown, as is a hypothetical signal with mass 
mvLQ'i = 320 GeV/c 2 , normalized to the data sample size. 



systematic uncertainty, as measured using methods de- 
scribed elsewhere [15j . is 3.0% of the selection efficiency 
and is primarily due to uncertainties in the response of 
the hadronic calorimeter to charged tracks and uncertain- 
ties in the track reconstruction efficiency. Uncertainty 
in simulating the effect of the lepton isolation require- 



ment is evaluated by comparing the efficiency of this iso- 
lation requirement in data and simulation. This yields a 
3.0% uncertainty on the selection efficiency. Additional 
sources of systematic uncertainty on the signal efficiency 
are the choices of the PDFs and the jet energy scale, as 
well as smaller contributions from the acceptance crite- 
ria, the electron (muon) identification, and the frp sim- 
ulation. The total systematic uncertainty on signal effi- 
ciency ranges from 10.5% (10.4%) for the er^ (/ttt,) chan- 
nel, for a VLQ3 with mass myLQ3 = 160GeV/c 2 , down 
to 7.0% (7.4%) for m VLQ3 = 400GeV/c 2 . The systematic 
uncertainties on the theoretical prediction of the VLQ3 
pair production cross section arise from the choices of 
PDFs and renormalization scale Q 2 . These two contri- 
butions are combined in quadrature. 

A likelihood function is constructed using a Poisson 
probability distribution of the expected rates of signal 
plus backgrounds and observed number of events in each 
channel and signal region. For each VLQ3 mass, the 
expected signal rates are functions of the VLQ3 pair pro- 
duction cross section. The expected signal rates also 
include factors for the branching ratios of leptonic and 
hadronic tau decays, luminosity, and full selection effi- 
ciencies. Systematic uncertainties, including 6% due to 
the luminosity measurement (not included in Table JJ, 
are incorporated into the fit. We apply Gaussian proba- 
bilities for the uncertainties on the background estimates, 
and account for correlations among different sources of 
systematic uncertainties. To set cross section limits for 
each mass, we integrate the likelihood distribution over 
all parameters except the cross section, and then inte- 
grate as a function of cross section from zero up to the 
cross section where the integral reaches 0.95. 

The results are shown in Fig. [2l as a function of 
VLQ3 mass, along with the theoretical predictions. For a 
VLQ3 with Yang-Mills couplings, at 95% confidence level 
(C.L.), the upper limit on the cross section is a < 344 fb, 
assuming B(VLQ3 — > br) = 1, and the lower limit on 
the mass is mvLQ3 > 317 GeV/ c 2 . With theoretical un- 
certainties included on the predicted cross section, the 
results are a < 360 fb and mvLQ3 > 294 GeV/c 2 . For the 
minimal couplings model, the upper limit on the cross 
section is a < 493 fb and the lower limit on the mass 
is mvLQ3 > 251 GeV/ c 2 . With theoretical uncertainties 
included on the predicted cross section, the results are 
a < 610 fb and myLQs > 223GeV/c 2 . The mass limits 
are approximately 80 — 90 GeV/c 2 higher than those of 
previous comparable results d, Q . 

Using 322 pb -1 of luminosity at CDF 11, we have 
searched for VLQ3 pair production and subsequent decay 
to two tau leptons and two jets. We observe no excess 
of events beyond the expected SM processes and set the 
most stringent limits to date on the VLQ3 mass and pair 
production cross section in the context of two coupling 
scenarios. 
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FIG. 2: The 95% C.L. upper limits for VLQ3 pair production 
versus mass. Also shown are the theoretical predictions using 
the simulation described in the text, with bands for uncer- 
tainties due to the choices of PDFs and Q 2 . 
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